The influence of multiple ionization thresholds on harmonic generation: Ar 
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We apply time-dependent _R-matrix theory to investigate harmonic generation from ground state 
Ar + with M = at a wavelength of 390-nm. Contributions associated with the different 3s 2 3p 4 
ionization thresholds are assessed, including the interference between these. The dominant contri- 
bution originates from the second ionization threshold, 3s 2 3p 4 1 D. Changes to the harmonic yields 
arising from the higher 3s3p 5 thresholds are also assessed. We further confirm that Ar + has a higher 
harmonic yield than He for the same laser pulse, despite having a higher ionization threshold. 

PACS numbers: 32.80.Rm, 31.15.A-, 42.65.Ky 



I. INTRODUCTION 

Harmonic generation (HG) has been at the centre of 
atomic and molecular physics for over 20 years Q. Not 
only is it the means of producing the ultrashort laser 
pulses which drive attosecond physics Q , but it has been 
adopted as an important measurement tool for some of 
the most fundamental physical atomic and molecular pro- 
cesses. HG has facilitated imaging molecular dynamics 
obtaining detailed information about molecular or- 
bitals [H and has been employed to highlight the impor- 
tance of electron correlation and multielectron dynamics 
in a variety of systems on ultrashort timescales [5j-l7j . 

Generally, HG is described in terms of a three-step 
model Q : a laser field causes tunnel ionization of an elec- 
tron, which is subsequently accelerated in the field. As 
the electric field changes direction, the electron is driven 
back towards its parent ion, and can be recaptured, emit- 
ting a high-energy photon. A semiclassical model of this 
process has proven highly successful in explaining exper- 
imental phenomena, particularly in noble-gas atoms 

The three-step model describes the main physics lead- 
ing to HG, but the model cannot be applied directly to all 
atomic or molecular systems. It assumes that tunnel ion- 
ization leaves the residual ion in the ground state. This 
is appropriate for noble gas atoms, but may not be so 
for systems in which the residual ionic state lies close to 
other ionic states. Examples of such systems are primar- 
ily found in molecules such as Ng. It is less appreciated 
that they can also be found in atomic systems, eg. Ar + . 
For these systems, electron-emission channels associated 
with different ionic states contribute to HG, and interfer- 
ence between these channels can be of significance. Since 
HG occurs within a laser field cycle, this interference pro- 
vides information about ultrafast dynamics. To develop 
understanding of how this information can be extracted, 
accurate theoretical data is of great benefit. To obtain 
such data, it is imperative to apply theoretical methods 
capable of including multiple channels associated with 
different thresholds and the interactions between these 
channels: recent studies on the Cooper minimum in Ar 
have shown that the inclusion of multichannel effects can 
alter the harmonic yield by as much as two orders of 



TABLE I: Energies of the five ionization thresholds of Ar + 
included in the present work with respect to the Ar 2+ ground 
state, and compared to literature values [14 ]. 



Configuration Term 



Energy (Lit. 
eV 



Energy (TDRM) 
eV 



3s 2 3p 4 



3s3p 5 



3p 

1 D 

3 po 
1 po 



0.00 
1.67 
4.06 
14.10 
17.79 



0.00 
2.04 
3.85 
17.04 
24.83 



magnitude [lOj . 

We have recently developed capability within time- 
dependent i?-matrix theory (TDRM) [ll| to determine 
harmonic spectra for general multielectron atoms 0, [HJ . 
TDRM is well-suited to study ultrafast multielectron dy- 
namics, as demonstrated in a study of C + [l3j, and mul- 
tichannel interferences in HG, as evidenced in a study 
of resonant enhancement of the 5th harmonic in Ar [7j . 
In order to study the interplay of multiple channels as- 
sociated with low-lying thresholds, we apply TDRM to 
investigate HG from Ar + in a 4 x 10 14 Wcm -2 , 390-nm 
laser field. This intensity /wavelength regime can lead to 
significant non-perturbative changes in the atomic struc- 
ture. Therefore, perturbative methods are not suitable 
for addressing an ion in these fields, even though the 
interaction can be characterised as a multiphoton pro- 
cess. Given the small energy gap between the relevant 
thresholds, we would expect that the general multichan- 
nel effects seen for 390-nm would still be important at 
800-nm, even though the fine detail may differ. 

Ar + is an ion well-suited to investigations on multiple 
thresholds: the lowest three ionization thresholds, 3s 2 3p 4 
3 P, 1 D and X S, are separated from each other by about 2 
eV (Tab. []}. We can also investigate the effects of higher- 
lying thresholds by including the 3s3p 5 thresholds. Ar + 
thus allows the investigation of interference effects arising 
from the interplay between channels associated with well 
separated, as well as closely spaced, thresholds. 

The harmonic response of Ar + to intense laser light is 
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also of relevance to experiment. The highest harmonics 
(up to 250 eV) generated by irradiation of Ar by intense 
laser light have been assigned to the response of Ar + 
[TBI Il6j . A full picture of HG at high intensities should 
therefore include ionized species. Recent photoionization 
experiments and calculations also provide detailed infor- 
mation on resonances that may affect HG in Ar + [I?} • 

II. THEORY 

The TDRM approach is a non-perturbative, ah ini- 
tio approach to describe multielectron atoms in short, 
intense li ght pulses. Full details of the method can be 
found in (ill ]. By propagating the solution of the time- 
dependent Schrodinger equation on a discrete time mesh 
of step size At, we can express its solution, +1 , at time 
t q +i in terms of the solution at the previous time step t q : 

(H m -E)* tq+1 =-(H m +E)* tq . (1) 

Here, H m is the Hamiltonian at the midpoint of the time 
interval, and it contains both the non-relativistic field- 
free atomic Hamiltonian and the laser interaction term. 
The laser field is assumed to be linearly polarized and 
spatially homogeneous and is described by the dipole ap- 
proximation in the length form jl8(. The imaginary en- 
ergy, E, is defined as 2i/At. 

TDRM theory makes use of the i?-matrix partition of 
configuration space. Within an inner region close to the 
nucleus, full account is taken of all electron-electron in- 
teractions. Outside of this region, exchange interactions 
between an ejected electron and those remaining close to 
the atomic core can be neglected, and the ejected elec- 
tron moves only in the long-range multipole potential of 
the residual ion and the laser field. 

We evaluate Eq. (JXJ) at the boundary of this inner 
region, ai n , as a matrix equation (Tlj : 

F(a in ) = R(a in )F(a in ) + T(a in ), (2) 

in which the wavefunction, F, at the boundary is ex- 
pressed in terms of its derivative, F, and an inhomoge- 
neous vector, T, which arises from the right hand side 
of Eq. (JT]) . The matrix R connects the inner and outer 
region wavefunction at the boundary, a in . 

Given an inner region wavefunction, R and T are eval- 
uated at the boundary, and propagated outwards in space 
to an outer region limit at which it can be assumed that 
the wavefunction has vanished. There, the wavefunction, 
F, can be set to zero and propagated inwards to a[ n . 
Once F is determined at every boundary point, the full 
wavefunction can be extracted from the i?-matrix equa- 
tions. We can then iterate the procedure using Eq. ((T|). 

The light radiated by an oscillating dipole is propor- 
tional to its acceleration [l9|, [2(| ■ It is also possible to 
express the harmonic spectrum in terms of the dipole ve- 
locity or dipole operators (2l| - [23 |. For He, the different 
expressions have been shown to be both self-consistent 



within TDRM, and consistent with the spectrum ob- 
tained by the HELIUM approach [HHl]. We show only 
dipole length spectra here. We have verified that they are 
consistent with those based on the dipole velocity form. 

The harmonic response of a single atom can be evalu- 
ated via the expectation value of the dipole operator: 

d(t)oc<*(t)M*(t)>, (3) 

where z is the total position operator along the laser 
polarization axis. The harmonic spectrum is then pro- 
portional to the square of the modulus of the Fourier 
transform of d(t), d(u;)| 2 . We note that this is a non- 
relativistic approximation to the HG process. Inclusion 
of relativistic effects, such as spin-orbit coupling, would 
lead to the population of M — 1 levels, and this is the 
subject of ongoing investigation. 



III. COMPUTATIONAL SETUP 

In i?-matrix theory, Ar + is described as a state of 
Ar 2+ plus an additional electron. To describe Ar 2+ , we 
use a set of Hartree-Fock {Is, 2s, 2p, 3s,3jp| orbitals, ob- 
tained for the Ar 2+ ground state [2g|. We obtain the 
3s 2 3p 4 and 3s3p 5 eigenstates of Ar 2+ from configuration- 
interaction calculations comprising 3s 2 3p 4 , 3s3p 5 and 
3p 6 . No pseudo-orbitals are included, since such orbitals 
may lead to spectra influenced by unphysical resonances. 
A consequence of this basis-set restriction is that the en- 
ergies of the Ar 2+ thresholds differ by 0.2 and 0.4 eV 
from experiment for the 3s 2 3p 4 states and 3 and 7 eV for 
the 3s3p 5 states (Tab. IJ). To assess interference effects, 
we employ several Ar + models. The full Ar + model con- 
tains all five Ar 2+ thresholds. In the three-state model, 
only the 3s 2 3p 4 states of Ar 2+ are retained. We also use 
Ar + models in which only a single state or a pair of states 
from the 3s 2 3p 4 configuration is retained. In these latter 
models, both 3s3p 5 states are always included. 

The inner region has a radius of 15 a.u. which suffices 
to contain the residual Ar 2+ ion. The outer region radius 
is 600 a.u. The set of continuum orbitals is built from 60 
£?-splines for each angular momentum of the continuum 
electron. The Ar + basis contains all allowed combina- 
tions of Ar 2+ states and the set of continuum orbitals up 
to a total angular momentum L max = 19. Convergence 
testing was carried out up to L max = 23. The outer 
region is divided into sectors of 2 a.u. each containing 
35 .B-splines of order 9 per channel. The time step in 
the wavefunction propagation is 0.1 a.u. We use 390-nm 
laser pulses, consisting of a 3 cycle sin 2 ramp-on followed 
by 2 cycles at peak intensity, 4 x 10 14 W/cm 2 , and a 3 
cycle sin 2 ramp-off. The initial state is the Ar + ground 
state with total magnetic quantum number M — 0, which 
would be the dominant non-relativistic Ar + ground state 
level following strong-field ionization of Ar at 390-nm. 
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FIG. 1: (Color online) Harmonic spectrum of Ar + produced 
by an 8-cycle 4 x 10 14 Won" 2 laser pulse at 390-nm, as ob- 
tained by the three-state model (dotted, red line) and by the 
five-state model (solid, blue line). 



IV. RESULTS 

Figure Q] shows the harmonic spectra produced from 
the full and the three-state model of Ar + . The three- 
step model [§[ suggests a harmonic cut-off of 46 eV. Even 
though the three-step model may not apply at 390-nm, 
the harmonic spectra are not inconsistent with this value. 

As shown in Fig. [T] inclusion of the 3s3p 5 thresh- 
olds has a noticeable effect on the yields for the 11th to 
the 15th harmonic, which are reduced by up to an order 
of magnitude. These harmonics coincide with the range 
of energies associated with the Rydberg series converg- 
ing onto the 3s3p 5 thresholds. We can thus ascribe the 
differences between the two spectra as arising from the 
reaction of multiple electrons to the laser field. Rydberg 
series converging onto 3s3p 5 effectively describe the ex- 
citation of a 3s electron, while the main contribution to 
the full harmonic spectrum involves the emission of a 3p 
electron. This competition is similar to the one observed 
for HG in Ar [7]. However, the present study employs 
a longer wavelength and involves a five-photon gap be- 
tween the 3s 2 3p 4 and 3s3p 5 thresholds. 

The main focus of the present study is interference 
between channels associated with the three low-lying 
3s 2 3p 4 thresholds. We have therefore calculated har- 
monic spectra using Ar + models in which only selected 
3s 2 3p 4 thresholds are included. 

Figure [2] compares the harmonic spectra obtained by 
summing contributions of the individual 3s 2 3p D and 
1 S models, the individual 3s 2 2p 4 3 P, 1 D and 1 S models 
and the spectrum from the full model. The individual 
3s 2 3p 4 models retain both 3s3p 5 thresholds. The spec- 
trum is largely dominated by the 1 D model, while the 3 P 
model does not contribute significantly, especially in the 
cutoff region, as demonstrated by the small difference 
between the two summed spectra. The largest differ- 
ence is seen for the 7th harmonic for which the 3 P model 
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FIG. 2: (Color online) Ar + harmonic spectrum produced by 
an 8-cycle, 4x 10 14 Wcm -2 laser pulse at 390-nm as calculated 
by summing contributions of the individual 3s 2 3p 4 1 D and 1 S 
models (dotted, black line), the individual 3s 2 2p 4 3 P, 1 D and 
S models (dashed red line) and from the full model (solid, 
blue line). All calculations include both 3s3p J thresholds. 

contributes 41%. In terms of the three-step mechanism, 
this suggests that tunnel ionization leaving Ar 2+ in the 
ground state does not significantly contribute to HG. In 
the three-step model for Ar + , the first step should con- 
sider tunnel ionization leaving Ar 2+ in an excited state. 

Figure [2] shows that interference between channels as- 
sociated with different 3s 2 3p 4 thresholds must be ac- 
counted for. The comparison between the summed spec- 
tra and the full model shows differences for the 3rd -7th 
harmonic and the 11th harmonic by as much as an or- 
der of magnitude, noticeable shifts in the energy of the 
13th and 15th harmonics and reductions by two orders of 
magnitude for the 21st and 23rd harmonics. The accu- 
rate determination of Ar + harmonic yields thus requires 
calculations including all 3s 2 3p 4 channels simultaneously. 

To assess the interferences, we have performed calcu- 
lations in which pairs of 3s 2 3p 4 thresholds of Ar 2+ arc 
retained. Of particular interest are the spectra in which 
we retain (a) 3s 2 3p 4 3 P and 1 D and (b) 3s 2 3p 4 3 P and 
1 S. As shown in Fig. [31 the ( 3 P, 1 £>) model provides a 
spectrum with harmonic peaks within 15% of the full 
model up to the 15th harmonic, apart from the 11th 
with a difference of 30%; the 17th and 19th peaks dif- 
fer by 25%. Channels associated with the 3 P and the 
1 D thresholds are hence the most important channels 
for HG. The improvement over the individual-state mod- 
els (Fig. [2]) shows that channels associated with the 3 P 
threshold are important to HG beyond the 3 P threshold, 
even though their direct contribution to these harmonics 
is relatively minor. 

Addition of the harmonic spectra obtained from the 
( 3 P, 1 Z?) and the ( 3 P, 1 S') model shows poor agreement 
for harmonics below the Ar + ionization threshold as the 
Rydberg series leading up to the 3 P threshold contributes 
twice. However, beyond this threshold, the agreement 
with the full model is very good, with harmonic peaks 
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FIG. 3: (Color online) Ar + harmonic spectrum produced by 
an 8-cycle 4 x 10 14 Wcm -2 laser pulse at 390-nm as obtained 
by the ( 3 P, 1 D) model (dotted, black line), by summing the 
spectra obtained by the ( 3 P, 1 D) model and the ( 3 P, 1 S) 
model (dashed, red line) and by the full model (solid, blue 
line) . 

differing by up to only 6% for the 9th to the 19th, apart 
from differences of 15% and 19% for the 11th and 17th 
harmonics. This comparison again indicates that HG 
due to channels associated with the 3 P threshold pri- 
marily occurs below the ionization threshold. Above this 
threshold, HG from these channels is negligible. 

The importance of the 3 P threshold for an accurate 
harmonic spectrum can be understood by considering the 
atomic structure. The first step within the three-step 
model is emission of a single electron. For an accurate 
emission rate, the atomic structure up to the threshold 
must be described well. Inclusion of the 3 P threshold 
is necessary to account for the Rydberg series leading 
up to the first ionization threshold. This series also af- 
fects the position of low-lying states such as 3s3p 6 and 
of low-lying members of Rydberg series leading up to 
other thresholds. Therefore, when HG originates from 
an excited threshold, atomic structure associated with 
the lower thresholds must still be accounted for. 

Finally, we compare the efficiency of HG in Ar + with 
that of other atoms. Ar + has a slightly higher ionization 
potential, 27.6 eV, than He, 24.6 eV. On the other hand, 
Ar + is a larger ion, and may therefore provide a greater 
harmonic response [l5| . We also consider the Ne + ion as 
it also is an ionized noble-gas atom, but with a substan- 
tially higher ionization potential, 40.96 eV. 

Figure 2] shows the harmonic spectrum of Ar 4 ", He and 
Nc + for the laser pulse used throughout. For He, we use 
the 6P-model described in our earlier work [l2j], whereas 
for Ne + we include the three 2s 2 2p 4 and the two 2s2p 5 
thresholds with the Ne 2+ states generated following a 
similar procedure as outlined above for Ar 2+ . Figure 0] 
shows expected behaviour for Ne + and its higher ioniza- 
tion potential, with a harmonic yield several orders of 
magnitude smaller than that of He and Ar + . However, 
the figure also shows that, apart from the 13th and 15th 



harmonics, the harmonic yield from Ar 4 " is consistently 
higher than that from He, despite its larger ionization 
potential. Thus Ar + is indeed an efficient ion for HG. 



V. CONCLUSION 

We have applied TDRM theory to determine harmonic 
yields for ground-state Ar + with M = and assess in 
detail the role of the various closely spaced ionization 
thresholds. The dominant contribution to the harmonic 
yield is associated with the first excited threshold instead 
of the lowest ionization threshold. This lowest threshold 
must still be accounted for, as it affects the atomic struc- 
ture leading up to excited thresholds. The 3s3p 5 thresh- 
olds affect the harmonic yield significantly less, but need 
to be taken into account for harmonics with energies close 
to these thresholds. Overall, the harmonic yield for Ar + 
is generally larger than the yield for He, even though Ar + 
has a larger ionization potential. 

In our discussion of the HG process we have frequently 
made reference to the three-step model, even though it 
may not be applicable this wavelength/intensity regime. 
However, the present calculation demonstrates that in- 
teraction between different channels is important for an 
accurate description of HG. Since the energy separation 
between the 3 P and 1 D thresholds is comparable to 800- 
nm photon energies we would expect these interactions 
to be important at these wavelengths as well. In these 
M = calculations, the highest harmonics (E > I p ) are 
strongly associated with the Ar 2+ 1 D channels. These 
channels would not be strongly disfavored at 800-nm be- 
cause of the small energy gap between the 3 P and 1 D 
thresholds. We would therefore expect the 1 D channels 
to remain important for HG in Ar + with M — at 800- 
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FIG. 4: (Color online) The harmonic spectrum produced by 
a 4 x 10 14 Wcm -2 390nm laser pulse interacting with Ar + 
(solid, red line), Ne + (dotted, black line) and He (dashed, 
blue line). The Ne + spectrum is orders of magnitude lower 
than the Ar + and He spectra, but the Ar + yield is slightly 
larger than the He, despite its higher ionization potential. 
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nm. Hence, application of the three-step model for Ar + 
with M = should account for multichannel interactions. 

The harmonic spectra obtained in the present study 
demonstrate a significant influence from interference as- 
sociated with different ionization thresholds. The use of 
multielectron codes, such as TDRM, is therefore essential 
for systems with excited thresholds just above the lowest 
ionization threshold. Atomic systems, such as Ar + , may 
be suitable for developing understanding of how mul- 
tichannel interactions affect HG. Further experimental 



studies of HG in Ar + would thus be very interesting. We 
have only addressed the non-relativistic case of M = 
here, but aim to extend the work to include the effect of 
M = 1 in our calculations. We also aim to extend these 
studies to longer wavelengths, which will enable more de- 
tailed comparisons with experiment. The recently devel- 
oped i?-matrix with time-dependence approach [13, 
will be of major benefit for this extension. 

ACB acknowledges support from DEL (NI). HWH is 
supported by EPSRC under grant number G/055416/1. 
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